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Abstract 


This  study  examines  the  ability  of  the  CTH  hydrocode  to  predict  die  effect  of  rod  nose-shape 
on  the  transition  from  rigid  body  to  eroding  rod  penetration  for  tungsten  alloy  long-rod 
penetrators  perforating  finite  aluminum  targets.  Two  rod  nose-shapes  and  two  target  alloys  were 
considered.  The  rod  nose-shapes  were  hemispherical  and  ogival,  and  the  target  alloys  were  7.62- 
cm-thick  5083  and  7039  aluminum.  Results  are  compared  to  an  experimental  study  that 
delineated  the  effect  of  nose-shape  on  the  threshold  velocity  at  which  tungsten  alloy  penetrators 
transition  from  rigid  body  to  eroding  rod  when  perforating  finite  aluminum  targets. 


11 


Acknowledgments 


The  author  would  like  to  thank  Dr.  Steven  B.  Segletes  and  Mr.  Kent  D.  Kimsey  for  their  helpful 
comments  and  suggestions  regarding  this  paper.  The  author  would  especially  like  to  thank  Dr.  Lee 
S.  Magness,  Jr.  for  discussions  regarding  his  experiments  and  assuring  the  acciuracy  of  the 
discussions  of  them.  Additional  thanks  go  to  Dr.  Stewart  A.  Silling  of  Sandia  National  Laboratories 
for  discussions  regarding  his  boundary  layer  interface  algorithm  and  for  providing  an  advanced  copy 
of  parts  of  his  user  manual. 


m 


« 


Intentionally  left  blank. 


IV 


Table  of  Contents 


Page 

Acknowledgments  . iii 

List  of  Figures  . vii 

List  of  Tables  .  ix 

1.  Introduction .  1 

2.  Problem  Setup  .  3 

3.  Results  and  Discussion .  6 

3.1  Residual  Shapes  for  Baseline  Simulations .  8 

3.2  Residual  Velocity  for  Baseline  Simulations  .  14 

3.3  BLINT  Model  Parameters .  17 

3.4  Failure  Model  Effects .  20 

4.  Conclusions .  22 

5.  References .  29 

Appendix  A:  Input  for  Hemi-Nose  KE-BL-CO-P  Sims  vs.  5083  Aluminum 

Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes . .  33 

Appendix  B:  Input  for  Ogival-Nose  KE-BL-CO-P  Sims  vs.  5083  Aluminum 

Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes  . .  41 

Appendix  C:  Input  for  Hemi-Nose  KE-BL-CO-P  Sims  ys.  7039  Aluminum 

Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes . .  49 

Appendix  D:  Input  for  Ogival-Nose  KE-BL-CO-P  Sims  vs.  7039  Aluminum 

Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes  . .  57 

Appendix  E:  Input  for  Hemi-Nose  KE-BL-CO-S  Sims  vs.  5083  Aluminum 

Tai^ets  .  65 

Appendix  F:  Input  for  Ogival-Nose  KE-BL-CO-S  Sims  vs.  5083  Aluminum 

Targets  .  73 


V 


Page 

Appendix  G:  Table  of  Simulation  Results  .  81 

Distribution  List .  85 

Report  Documentation  Page .  91 


VI 


List  of  Figures 


Figure  Page 

1.  Penetrator  Geometries .  3 

2.  Comparison  of  Hemi-Nose  Penetrator  KE  and  MV  Sims  Residual  Shapes  With 

Experiment  After  Perforating  5083  Aluminum  Targets  .  9 

3.  Comparison  of  Hemi-Nose  Penetrator  KE  and  MV  Sims  Residual  Shapes  With 

Experiment  After  Perforating  7039  Aluminum  Targets  .  10 

4.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  KE  Sims  After  Perforating 

5083  Aluminum  Targets  .  12 

5.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  MV  Sims  After  Perforating 

5083  Aluminum  Targets  .  13 

6.  Experimentally  Determined  Residual  Shapes  of  the  Ogival-Nose  Penetrator  After 

Perforating  7039  Aluminum  Targets .  14 

7.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  KE  Sims  After  Perforating 

7039  Aluminum  Targets  .  15 

8.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  MV  Sims  After  Perforating 

7039  Aluminum  Targets  . .  16 

9.  Residual  Velocity  Comparison  Between  Experiment,  KE,  and  MV  Sims  After 

Perforating  5083  Aluminum  Targets .  17 

10.  Residual  Velocity  Comparison  Between  Experiment,  KE,  and  MV  Sims  After 

Perforating  7039  Aluminum  Targets .  18 

1 1 .  Residual  Shapes  of  KE-NB-NC-P  Sims  After  Perforating  5083  Aluminum 

Targets .  19 

12.  Comparison  of  Residual  Velocity  for  KE-NB-NC-P  Sims  With  Experiment  and 

Baseline  KE  Sims  After  Perforating  5083  Aluminum  Targets .  20 

13.  Residual  Shapes  for  Hemi-Nose  Penetrator  KE-BL-NC-P  Sims  After 

Perforating  5083  Aluminum  Targets .  21 


vii 


Figure  Page 

14.  Residual  Shapes  for  Ogival-Nose  Penetrator  KE-BL-NC-P  Sims  After  Perforating 

5083  Aluminum  Targets  .  22 

15.  Comparison  of  Residual  Velocity  for  KE-BL-NC-P  Sims  With  Experiment  and 

Baseline  KE  Sims  After  Perforating  5083  Aluminum  Targets .  23 

16.  Residual  Shapes  for  Hemi-Nose  Penetrator  KE-BL-CO-S  Sims  After  Perforating 

5083  Aluminum  Targets  .  24 

1 7.  Residual  Shapes  for  Ogival-Nose  Penetrator  KE-BI^CO-S  Sims  After  Perforating 

5083  Aluminum  Targets  .  25 

1 8.  Comparison  of  Residual  Velocity  for  KE-BL-CO-S  Sims  With  Experiment  and 

Baseline  KE  Sims  After  Perforating  5083  Aluminum  Targets .  26 


List  of  Tables 


Table  Page 

1.  Equation-of-State  Parameters  .  4 

2.  Initial  Impact  Conditions  and  Ballistic  Test  Results .  7 


IX 


Intentionally  left  blank. 


1.  Introduction 


In  examining  the  convergence  characteristics  of  the  Eulerian  CTH  hydrocode  [1]  as  a  function 
of  spatial  resolution,  Zukas  [2]  found  that  the  code  could  not  accurately  predict  perforation  of  armor 
plate  by  a  hard  projectile  at  low  velocities  (less  than  1.5  km/s).  A  penetrator,  which,  in  experiments, 
perforated  a  finite  steel  target  with  significant  residual  length  and  velocity,  was  predicted  to  be 
unable  to  perforate  the  target.  Previously,  this  problem  had  been  modeled  successfully  using  an 
in-house  version  of  the  EPIC  Lagrangian  hydrocode.  Zukas  observed  that,  regardless  of  the  mixed 
cell  strength  formulation  used  (several  are  available  in  the  CTH  hydrocode),  high-strength  penetrator 
material  included  in  a  mixed  cell  was  modeled  as  being  significantly  softer — ^an  unrealistic  treatment 
that  caused  excessive  deformation  in  the  penetrator.  The  net  effect  was  that  the  CTH  hydrocode 
could  not  accurately  model  the  rigid-body  penetration  of  a  soft  target,  an  eroding  projectile 
penetrating  harder  targets  at  low  velocities  or  the  sliding  between  two  material  interfaces. 

A  new  boundary  layer  algorithm  for  sliding  interfaces  (BUNT)  was  recently  incorporated  into 
the  CTH  hydrocode  for  two-dimensional  problems  only  [3].  The  algorithm  relocates  the  slip  layer 
outside  of  mixed  cells  into  the  softer  material,  thus  allowing  hard  materials  to  penetrate  as  rigid 
bodies.  Good  correlation  with  experiments  has  been  obtained  using  the  BUNT  algorithm  by  Silling 
[4]  and  Kmetyk  and  Yarrington  [5].  Both  modeled  hard  paietrators  impacting  soft  targets  knowing, 
a  priori,  that  the  penetrators  would  remain  rigid. 

This  study  examined  the  ability  of  the  CTH  hydrocode  (August  1993  release)  to  predict  the 
impact  velocity  at  which  a  penetrator  would  transition  from  rigid-body  to  eroding-rod,  the  effect  of 
the  rod’s  nose  shape  on  this  transition  velocity  and  its  residual  velocity  and  shape  while  perforating 
finite  aluminum  targets.  The  perforation  of  soft  aluminum  targets  by  tungsten  alloy 
(95W-2.5Ni-l.0Fe-l.5Co,  cold  worked  by  swaging  to  a  21%  reduction  in  area)  long  rods  was 
modeled.  To  gauge  the  accuracy  of  the  CTH  hydrocode  with  the  BUNT  algorithm,  the  simulation 
results  were  compared  to  the  experimentally  determined  residual  penetrator  velocity  and  shape  [6]. 
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Two  rod  nose-shapes  and  two  target  alloys  were  considered.  The  rod  nose-shapes  were 
hemispherical  and  ogival,  and  the  target  alloys  were  7.62-cm  (3  in)-  thick  5083  and  7039  aluminum. 

A  subset  of  the  results  of  this  study  first  appeared  at  the  Symposium  on  Structures  Under 
Extreme  Loading  Conditions  as  part  of  the  1996  ASME  Pressure  Vessel  and  Piping  Conferaice  that 
took  place  in  Montreal,  Canada  from  21-26  July  1996  [7].  A  companion  paper  [8]  and  report  [9] 
also  exist,  which  provide  the  simulation  results  for  impacts  with  53.34-cm-thick  aluminum  targets 
where  results  are  compared  with  experimental  depth-of-penetration  (DOP)  tests.  The  companion 
report  provides  more  detail  than  the  companion  paper.  This  report  differs  from  the  original  paper 
by  including  additional  simulations  and  experiments,  a  more  detailed  discussion  of  the  BUNT 
model,  the  examination  of  effective  plastic  strains,  as  well  as  the  input  decks  used  for  the 
simulations. 

The  CTH  hydrocode  is  a  state-of-die-art,  second-order  accurate,  Eulerian  hydrocode  developed 
by  Sandia  National  Laboratories  which  is  capable  of  solving  complex  problems  in  shock  physics  in 
one,  two,  or  three  dimensions.  The  code  provides  several  constitutive  models,  including  an 
elastic-perfecdy  plastic  model  with  provisions  for  work  hardening  and  thermal  softening,  the 
Johnson-Cook  model  [10],  the  Zerrilli-Armstrong  model  [11],  the  Steinberg-Guinan-Lund 
model  [12, 13],  and  an  undocumented  power-law  model.  High-explosive  detonation  can  be  modeled 
using  the  programmed  bum  model,  the  Chapman-Jouguet  volume  bum  models,  or  the  history 
variable  reactive  bum  model  [14].  Several  equation-of-state  (EOS)  options  are  available,  including 
tabular  (i.e.,  SESAME),  analytical  (ANEOS),  Mie-Griineisen,  and  Jones-Wilkins-Lee  (JWL)  [15]. 
Material  failure  occurs  when  a  threshold  value  of  tensile  stress  or  hydrostatic  pressure  is  exceeded. 
In  addition,  the  Johnson-Cook  failure  model  [16]  is  also  available.  When  failure  occurs  in  a  cell, 
void  is  introduced  until  the  stress  state  of  the  cell  is  reduced  to  zero.  Recompression  is  permitted. 
To  reduce  the  diffusion  typically  encountered  in  Eulerian  simulations,  several  advanced  material 
interface  tracking  algorithms  are  provided,  including  the  high-resolution  interface  tracking  (HRIT) 
algorithm  (available  for  two-dimensional  simulations  only),  the  simple  line  interface  calculation 
(SUC)  algorithm  [17],  and  the  Sandia-modified  Young’s  reconstraction  algorithm  (SMYRA)  [18]. 
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2.  Problem  Setup 


The  two  geometries  for  the  tungsten  alloy  penetrators  are  shown  in  Figure  1.  Both  of  the 
penetrators  have  a  length  of  10.1346  cm  (3.99  in)  and  a  diameter  of  0.67564  cm  (0.266  in).  Due  to 
their  different  nose  shapes,  the  masses  of  the  penetrators  differ  slightly.  The  mass  of  the  hemi-nose 
penetrator  is  approximately  65  g,  and  that  of  the  ogival-nose  penetrator  is  approximately  63  g. 


10.1346 


0.6756 


1 

1 


All  Dimensions  Are  In  Centimeters 

Figure  1.  Penetrator  Geometries. 

Three  different  constitutive  models  were  used  in  the  simulations  to  model  the  deviatoric  response 
of  the  materials.  The  choice  of  the  constitutive  model  used  for  a  material  was  governed  by  the 
availability  of  material  data.  Material  data  were  not  available  for  the  95W-2.5Ni-l.0Fe-l.5Co,  21% 
swaged  tungsten  aUoy  penetrators  used  in  the  experiments.  Therefore,  the  alloy  was  approximated 
using  95W-3.5Ni-l.5Fe  tungsten  aUoy  data  for  the  Steinberg-Guinan-Lund  strain-rate-independent 
model  reported  in  Steinbeig  [19].  This  tungsten  alloy  has  die  same  percentage  of  tungsten  and  samR 
approximate  density  as  the  95W-2.5Ni-l.0Fe-l.5Co,  21%  swaged  alloy.  For  the  7039  aluminum 
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target,  the  Johnson-Cook  constitutive  model  was  used  with  the  parameters  reported  in  Johnson  and 
Cook  [10].  For  the  5083  aluminum  target,  a  power-law  constitutive  model  was  used  with  the 
parameters  reported  in  Silling  [3]  and  originally  reported  in  Forrestal  et  al.  [20]. 

The  Mie-Griineisen  EOS  was  used  for  all  materials.  EOS  data  were  obtained  from  a  data  file 
provided  with  the  CTH  hydrocode.  The  EOS  parameters  for  5083  aluminum,  7039  aluminum,  and 
95%  tungsten  content  tungsten  alloy  were  not  available.  Therefore,  they  were  approximated  using 
parameters  for  6061  aluminum,  7075  aluminum,  and  90W-7Ni-3Fe  tungsten  alloy,  respectively.  The 
initial  densities  of  the  6061  and  the  7075  aluminum  alloys  were  changed  to  reflect  those  for  5083 
and  7039  aluminum,  as  reported  in  the  Metals  Handbook  Desk  Edition  [21].  The  initial  density  of 
the  90W-7Ni-3Fe  alloy  was  changed  to  reflect  the  initial  density  of  the  95W-3.5Ni-l.5Fe  alloy 
reported  in  Steinberg  [19].  The  EOS  parameters  used  for  the  materials  are  listed  in  Table  1. 


Table  1.  Equation-of-State  Parameters 


Material 

Density 

Sound  Speed 

Slope  Us-Up 

Gruneisen 

Specific  Heat 

Po 

Co 

Parameter 

(g/cm^) 

(km/s) 

(s) 

(r„) 

(erg/g/eV) 

WAUoy 

18.16 

4.03 

1.237 

1.67 

1.66el0 

5083  Al 

2.66 

5.34 

1.40 

1.97 

1.07ell 

7039  Al 

2.77 

5.20 

1.36 

2.20 

1.07ell 

Failure  in  most  of  the  simulations  was  modeled  using  a  threshold-pressure  criterion.  The  tensile 
pressure  at  which  the  tungsten  alloy,  the  5083  aluminum,  and  the  7039  aluminum  were  assumed  to 
fail  was  3.5, 0.45,  and  0.50  GPa,  respectively.  Additional  simulations  used  a  strain-based  failure 
criterion  that  is  described  later  in  this  report. 

All  simulations  used  a  two-dimensional  cylindrical  coordinate  mesh  consisting  of  85  x  832  cells. 
The  mesh  in  the  radial  direction  starts  at  the  axis  of  symmetry  with  a  constant  cell  size  of 
0.0422275  cm  out  to  a  radius  of  1.6891  cm.  Thereafter,  cell  dimensions  expand  by  5%  increments 
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out  to  the  outer  radius  of  the  target.  This  mesh  provides  eight  cells  across  the  radius  of  the 
penetrator.  The  mesh  in  the  axial  direction  mesh  has  a  constant  cell  size  of  0.0422275  cm.  Thus, 
cells  in  the  penetrator-target  interaction  region  have  a  one-to-one  aspect  ratio. 


Parameters  for  the  BUNT  model  were  chosen  to  be  similar  to  those  reported  in  Kmetyk  and 
Yarrington  [5].  Thus,  the  boundary-layer  distance  (Wh)  and  the  slip-layer  distance  (Wj,)  were  chosen 
to  be  twice  the  zone  size  of  cells  in  the  penetrator-target  interaction  region.  The  boundary-layer 
distance  defines  which  cells  will  be  included  in  the  boundary  layer,  ff  the  cell  center  of  a  cell  is 
located  away  from  a  cell  whose  center  is  included  in  the  interface  layer,  it  is  considered  to  be  part 
of  the  boundary  layer.  Materials  defined  as  “hard”  make  up  the  hard  boundary  layer  and  materials 
defined  as  “soft”  make  up  the  soft  boundary  layer.  The  interface  layer,  which  is  about  two  cell 
widths  thick,  contains  all  cells  whose  hard  and  soft  material  volume  vector  gradient  magnitudes  are 
both  greater  than  or  equal  to  0.1.  The  slip-layer  distance  defines  which  cells  will  be  included  in  the 
slip  layer.  If  the  cell  center  of  a  cell  located  in  the  soft  boundary  layer  is  Wji  from  a  cell  whose  center 
is  included  in  the  interface  layer,  it  is  considered  part  of  the  slip  layer.  Cells  located  in  the  slip  layer 
have  their  flow  stresses  set  to  zero,  allowing  sliding  to  occur  in  these  cells.  An  option  to 
automatically  increase  the  yield  strength  of  the  penetrator  material  by  a  factor  equal  to 
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(where  r^  is  the  outer  radius  of  the  penetrator)  was  used.  This  ratio  represents  the  cross-sectional  area 
of  the  penetrator  plus  the  boundary-layer  distance  over  the  original  penetrator  cross-sectional  area. 
The  option  was  used  because  numerical  noise  can  cause  shear  stresses  close  to  the  yield  stress  to 
exceed  the  yield  stress,  causing  premature  irreversible  deformation  of  the  penetrator.  An  additional 
option  allows  for  the  inclusion  of  friction;  however,  ftiction  between  the  target  and  penetrator  was 
not  accounted  for  in  this  study.  Kmetyk  and  Yarrington  [5]  showed  that  the  BUNT  model  tended 
to  oveipredict  penetration  in  deep  penetration  problems  unless  friction  was  included. 


The  CTH  hydrocode  (August  1993  version)  cannot  convect  velocity  in  a  manner  such  that  both 
momentum  (MV)  and  kinetic  energy  (KE)  are  both  conserved  exactly.  The  default  option  allows 


5 


conservation  of  KE  such  that  total  energy  is  conserved  during  the  convection  phase  of  a 
computational  cycle;  however,  MV  is  not  conserved.  A  second  option  convects  velocity  such  that 
MV  is  conserved  during  the  convection  phase  of  a  computational  cycle  and  any  KE  discrepancies 
are  discarded.  Simulations  were  run  for  both  of  these  convection  options.  A  final  option  conserves 
both  MV  and  total  energy  during  the  convection  phase  of  a  computational  cycle  by  depositing  the 
KE  discrepancy  into  internal  energy.  This  can  have  the  effect  of  artificially  heating  a  material  [22], 
and  therefore,  this  option  was  not  used.  (Note:  With  the  March  1995  release  of  the  code,  a  half 
index  shifted  momentum  scheme  was  introduced  as  the  only  convection  option,  thus  the  choice  of 
convection  options  discussed  previously  are  no  longer  available  in  code  versions  later  than  August 
1993.) 

Complete  listings  of  the  CTH  input  decks  used  for  the  simulations  are  given  in  Appendices  A-F. 

the  only  difference  in  the  input  decks  was  the  penetrator  striking  velocity,  conservation  method, 
or  whether  or  not  the  BLINT  model  with  or  without  the  strength  correction  factor  was  used,  then 
those  decks  are  not  listed.  Notes  in  the  input  decks  describe  the  required  changes  needed  for  the 
input  decks  not  listed. 

3.  Results  and  Discussion 

Initial  impact  conditions  and  ballistic  test  results  [6]  from  the  experiments  are  provided  in 
Table  2  and  simulation  results  are  provided  in  Appendix  G.  The  total  penetrator  yaw  in  the 
experiments  was  small  (in  most  cases  less  than  1  °).  However,  these  values  still  exceed  the  critical 
yaw  as  defined  in  Bjerke  et  al.  [23, 24]  since  the  penetration  channel  diameter  is  about  the  same  as 
the  penetrator  shank  diameter  for  rigid  body  penetration  of  soft  targets  (e.g.,  see  Fonestal  et  al.  [25]). 
Any  effects  of  yaw  were  not  treated  in  the  simulations. 

In  the  following  discussion,  simulations  will  be  designated  by  the  options  used.  The  designation 
will  be  of  the  form  XX-XX-XX-X,  where  the  first  set  of  X’s  represents  whether  or  not  KE  or  MV 
was  conserved  during  the  convection  phase  of  a  computational  cycle,  the  second  set  of  X’s 
represents  whether  the  BUNT  model  was  used  (BL)  or  not  used  (NB),  the  third  set  of  X’s  represents 
whether  the  yield  strength  correction  factor  was  used  with  the  BUNT  model  (CO)  or  not  used  (NC), 
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Table  2.  Initial  Impact  Conditions  and  Ballistic  Test  Results 


Shot  No. 

Total 

Yaw 

(0) 

Striking 

Velocity 

(m/s) 

Original 

Mass 

(g) 

Residual 

Velocity 

(m/s) 

Comments 

1  Hemi-Nose  Penetrators  vs.  7.62-cm-Thick  7039  Aluminum  || 

4326 

0.56 

1,198 

65.2 

1,020 

Eroded 

4327 

0.71 

1,038 

65.2 

964 

Rigid? 

4328 

0.71 

1,093 

65.5 

961 

Bulge,  Fractured 

1  Hemi-Nose  Penetrators  vs.  7.62-cm-Thick  5083  Aluminum  || 

4329 

0.25 

1,283 

65.7 

1,168 

Eroded 

4330 

0.35 

1,201 

65.7 

1,109 

Slight  Bulge 

4331 

0.25 

1,147 

65.8 

1,069 

Rigid,  Fractured 

II  Ogival-Nose  Penetrators  vs.  7.62-cm-Thick  5083  Aluminum  || 

4332 

0.25 

1,286 

63.2 

1,245 

Rigid 

4333 

0.56 

1,399 

63.5 

1,352 

Rigid 

4334 

0.35 

1,534 

63.5 

1,493 

Rigid 

4335 

0.00 

1,600 

63.6 

1,562 

Rigid 

1  Ogival-Nose  Penetrators  vs.  7.62-cm-Thick  7039  Aluminum  || 

4336 

1.03 

1,474 

63.3 

1,414 

Rigid,  Fractured 

4337 

0.56 

1,595 

63.5 

1,528 

Rigid,  Fractured 

4451 

0.00 

1,755 

-63.5 

1,677 

Slight  Bulge 

4452 

0.79 

1,768 

-63.5 

1,652 

Eroded,  Bent 

and  the  final  X  represents  whether  the  threshold  pressure  failure  model  was  used  (P)  or  the  strain- 
based  failure  model  was  used  (S).  For  example,  die  designation  KE-BI^CO-P  Sim  would  mean  that 
this  simulation  conserved  KE  during  the  convection  phase  of  a  computational  cycle,  that  the  BLINT 
model  was  used  with  the  strength  correction  factor,  and  that  the  threshold  pressure  was  used  to 
model  failure.  The  designations  KE-Bl^CO-P  and  MV-BI^CO-P  represent  the  baseline  simulations 
and  will  for  simplicity  also  be  called  KE  Sim(s)  or  MV  Sim(s),  respectively. 
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3.1  Residual  Shapes  for  Baseline  Simulations.  Figure  2  compares  the  hemi-nose  penetrator 
shapes  predicted  by  the  KE  and  MV  Sims  with  tracings  of  residual  penetrator  shapes  obtained  from 
radiographs  of  the  experiments  against  the  5083  aluminum  targets.  While  the  simulation  results  are 
all  to  the  same  scale  and  can  be  directly  compared  to  their  initial  geometry  (Figure  2a),  the 
experimental  penetrator  shapes  may  not  be  at  the  same  scale.  The  effective  plastic  strains  in  the 
penetrator  for  the  simulations  are  also  shown  with  strains  less  than  5%  not  being  plotted.  Visible 
deformation  of  the  penetrator  occurs  when  plastic  strains  between  45  and  55%  (represented  by  green) 
appear  on  axis  of  the  penetrator.  The  tracings  show  that  the  experimental  hemi-nose  penetrator 
remained  rigid  (although  its  tail  fractured)  at  a  striking  velocity  of  1,147  m/s  (Figure  2b).  At  a 
striking  velocity  of  1,201  m/s,  the  onset  of  plastic  deformation  of  the  nose  was  observed  (Figure  2c), 
and  at  a  striking  velocity  of  1,283  m/s,  the  penetrator  was  significantly  eroded  (Figure  2d).  The 
condition  of  the  hemi-nose  penetrator  predicted  by  the  KE  Sim  of  the  1 , 1 47-m/s  test  was  also  rigid, 
although  the  fracture  of  the  tail  observed  experimentally  did  not  occur  (Figure  2b).  Erosion  of  the 
penetrator  is  evident  in  the  KE  Sim  for  the  striking  velocity  of  1,201  m/s  (Figure  2c).  At  a  striking 
velocity  of  1 ,283  m/s,  the  residual  penetrator  predicted  by  the  KE  Sim  shows  about  the  same  amount 
of  nose  deformation  and  rod  length  as  in  the  experiment  (Figure  2d).  For  the  MV  Sims,  the 
hemi-nose  penetrator  also  remains  rigid  at  a  striking  velocity  of  1,147  m/s  (Figure  2b).  At  striking 
velocities  of  1,201  m/s  and  1,283  m/s  (Figures  2c  and  2d),  the  MV  Sims  displayed  much  less  plastic 
deformation  and  erosion  of  the  rod  nose  than  the  experiments  or  the  KE  Sims. 

The  final  hemi-nose  penetrator  shapes  from  both  the  experiments  and  the  baseline  simulations 
for  the  7039  aluminum  targets  are  shown  in  Figure  3.  The  apparent  yaw  seen  in  the  figure  for  the 
experiments  is  not  representative  of  any  yaw  the  penetrator  may  have  experienced  in  the 
experiments.  Experimentally,  the  hemi-nose  penetrators  fractured  at  the  two  lower  velocities  tested 
(Figures  3b  and  3c).  At  a  striking  velocity  of  1,038  m/s  it  is  unclear  whether  the  penetrator  in  the 
experiment  remained  rigid  (Figure  3b).  Erosion,  plastic  deformation,  and  bending  of  the  penetrator 
are  evident  at  a  striking  velocity  of  1,093  m/s  (Figure  3c).  At  a  striking  velocity  of  1,198  m/s,  the 
eroded  penetrator  shows  a  sizable  bulge  at  the  nose  as  well  as  yaw  into  or  out  of  the  page 
(Figure  3d).  The  penetrator  in  the  KE  Sims  remained  essentially  rigid  at  a  striking  velocity  of 
1,038  m/s  (Figure  3b).  Plastic  deformation  occurred  at  a  striking  velocity  of  1,093  m/s  (as  seen  by 
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Figure  2.  Comparison  of  Hemi-Nose  Penetrator  KE  and  MV  Sims  Residual  Shapes  With 
Experiment  After  Perforating  5083  Aluminum  Targets. 
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(a)  Initial  Geometry 
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Figure  3.  Comparison  of  Hemi-Nose  Penetrator  KE  and  MV  Sims  Residual  Shapes  With 
Experiment  After  Perforating  7039  Aluminum  Targets. 
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the  green  color  in  the  nose  in  Figure  3c).  A  large  bulge  and  erosion  in  the  penetrator  are  evident  at 
a  striking  velocity  of  1,198  m/s  (Figure  3d).  For  the  MV  Sims,  the  penetrator  also  remains  rigid  at 
a  striking  velocity  of  1,038  m/s  (Figure  3b).  The  predicted  length  and  shape  of  the  penetrator  at  a 
striking  velocity  of  1,093  m/s  were  about  the  same  as  they  had  been  at  a  striking  velocity  of 
1,038  m/s  (Figure  3c).  The  penetrator  is  clearly  deformed  at  a  striking  velocity  of  1,198  m/s  as  seen 
from  plastic  strains  in  the  green  range  at  the  axis  in  the  nose  of  the  penetrator  (Figure  3d). 

Because  the  KE  Sims  and  hemi-nose  options  seem  to  predict  the  onset  of  visible  deformation  and 
predict  rod  shape  better  then  the  MV  Sims,  most  additional  simulations  were  all  modeled  with  the 
KE  convection  conservation  option.  Results  for  the  ogival-nose  penetrator  MV  Sims  will  be 
presented  only  at  striking  velocities  for  which  experimental  data  were  available  at  the  time  of  the 
original  paper  [8];  therefore,  no  additional  MV  Sims  will  be  presented.  Predictions  of  the  ogival- 
nose  penetrator’ s  threshold  velocity  were  completed  in  advance  of  the  experiments.  Because  it  was 
felt  that  the  threshold  velocity  of  the  ogival-nose  penetrator  may  exceed  the  experimental  gun 
system’s  maximum  launch  velocity  and  because  of  funding  constraints,  experiments  were  not 
completed  for  the  5083  aluminum  target. 

Figures  4  and  5  show  the  residual  shapes  and  plastic  strains  for  the  ogival-nose  penetrator  after 
perforating  5083  aluminum  targets  for  the  KE  and  MV  Sims,  respectively.  The  experimental 
penetrator  shapes  are  not  shown  because  they  remain  rigid  at  all  striking  velocities  tested  (see 
Table  2).  It  can  be  seen  from  Figure  4  that  the  amount  of  plastic  stain  experienced  by  the  penetrators 
increased  with  striking  velocity.  Not  until  the  penetrator  experienced  plastic  strain  between  46  and 
55%  (beginning  of  the  green  range)  on  the  penetrator  axis  can  one  see  significant  difference  in  the 
residual  shape  of  the  penetrator  (Figure  4h).  From  Figure  4,  it  is  clearly  evident  that  visible  plastic 
deformation  occurs  at  striking  velocities  between  1,900  and  2,000  m/s  for  the  KE  Sims.  At  a  striking 
velocity  of  2,100  m/s,  the  KE  Sim  ogival-nose  penetrator  is  clearly  eroded  (Figure  4i).  The  MV  Sim 
ogival-nose  penetrators  of  Figure  5  seem  to  experience  larger  plastic  strain  at  the  same  corresponding 
striking  velocity  than  do  their  KE  Sim  counterparts  (Figure  4). 
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Figure  4.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  KE  Sims  After  Perforating  5083 
Aluminum  Targets. 
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Figure  5.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  MV  Sims  After  Perforating  5083 
Aluminum  Targets. 

Figures  6,  7,  and  8  show  the  residual  shapes  for  the  ogival-nose  penetrator  impacting  7039 
aluminum  targets  from  the  experiments,  KE  Sims  and  MV  Sims,  respectively.  For  the  experiments, 
only  two  shapes  are  shown  at  the  striking  velocities  between  which  the  transition  from  rigid-body 
penetration  to  eroding  rod  occurred.  These  two  experiments  were  conducted  after  the  simulations 
were  completed.  Residual  shapes  of  lower  velocity  ogival-nose  penetrator  experiments  are  not 
shown  as  they  all  remain  rigid.  In  the  experiments,  plastic  deformation  and  bending  of  the  penetrator 


(a)  1,755  m/s  Impact  Velocify 


(b)  1,768  m/s  Impact  Velocity 

Figure  6.  Experimentally  Determined  Residual  Shapes  of  the  Ogival-Nose  Penetrator  After 
Perforating  7039  Aluminum  Targets. 

are  evident  for  a  striking  velocity  of  1,755  (Figure  6a).  With  an  increase  in  striking  velocity  of 
13  m/s,  the  penetrator  is  clearly  eroded  (Figure  6b),  suggesting  that  the  transition  velocity  lies 
between  a  striking  velocity  of  1,755  and  1,768  m/s.  The  KE  Sims  predicted  that  the  onset  of  plastic 
deformation  occurred  between  a  striking  velocity  of  1,800  and  1,900  m/s  (Figures  7d  and  7e)  and 
that  erosion  occurred  between  a  striking  velocity  of  1,900  and  2000  m/s.  Again  the  MV  Sims  seem 
to  show  much  larger  plastic  strains  at  the  corresponding  striking  velocity  then  do  the  KE  Sims 
(Figures  7  and  8). 

3.2  Residual  Velocity  for  Baseline  Simulations.  Figure  9  compares  the  predicted  residual 
velocity  of  the  penetrators  to  the  experimentally  determined  residual  velocity  for  the  5083  aluminum 
targets.  The  experimental  results  are  represented  with  solid  symbols,  the  KE  Sims  are  represented 
with  hollow  symbols,  and  the  MV  Sims  are  represented  with  half-filled  symbols.  In  addition,  the 
hemi-nose  penetrators  are  represented  by  circles  and  the  ogival-nose  penetrators  are  represented  by 
squares.  In  all  cases,  the  predicted  residual  velocities  are  less  than  the  experimentally  determined 
ones.  Not  much  difference  is  seen  in  the  predicted  residual  velocities  of  the  ogival-nose  penetrator, 
though  the  KE  Sims  show  a  slight  improvement  over  the  MV  Sims.  Most  ogival-nose  predictions 
were  within  8.1%  of  the  experimentally  determined  residual  velocities  with  the  exception  being  the 
MV  Sim  at  a  striking  velocity  of  1 ,286  m/s,  which  differed  by  1 3. 1  %.  For  the  hemi-nose  penetrators. 
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Figure  7.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  KE  Sims  After  Perforating  7039 
Aluminum  Targets. 


the  predicted  residual  velocities  for  the  simulations  were  much  better  for  the  MV  Sims  than  for  the 
KE  Sims,  except  for  the  datum  at  the  lowest  striking  velocity.  In  general,  the  difference  between  the 
predicted  and  the  experimentally  determined  residual  velocities  for  the  hemi-nose  penetrators  was 
less  than  7.4%  for  KE  Sims  and  less  than  2.3%  for  MV  Sims,  except  for  the  damrn  at  the  lowest 
striking  velocity,  which  differed  by  9.0%. 
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Figure  8.  Residual  Shapes  of  the  Ogival-Nose  Penetrator  MV  Sims  After  Perforating  7039 
Aluminum  Targets. 


Figure  10  compares  the  experimental  residual  velocity  data  with  the  baseline  simulations’ 
predicted  residual  velocity  for  the  7039  aluminum  targets.  The  symbols  are  the  same  as  were  used 
in  Figure  10  for  the  5083  aluminum  target.  Again,  the  CTH  hydrocode’s  predicted  residual 
velocities  are  less  than  those  obtained  experimentally  for  all  cases  considered.  The  predicted 
residual  velocities  for  the  ogival-nose  penetrator  were  within  8.4%  of  those  obtained  experimentally, 
with  only  a  slight  difference  in  the  predictions  between  the  KE  and  MV  Sims.  The  fact  that  the 
hemi-nose  penetrator  experiments  show  almost  identical  residual  velocities  for  the  striking  velocities 
of  1,038  m/s  and  1,093  m/s  suggests  a  transition  from  rigid  body  penetrator  to  eroding/deforming 
penetration  between  these  two  velocities.  The  hemi-nose  KE  Sims  show  near  identical  residual 
velocities  for  striking  velocities  of  1,093  m/s  and  1,198  m/s  suggesting  that  the  transition  velocity 
lies  between  these  two  velocities.  For  the  MV  Sims,  it  is  not  clear  from  looking  at  the  residual 
velocities  whether  a  transition  from  rigid  body  to  eroding  rod  penetration  took  place.  The  predicted 
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Figure  9.  Residual  Velocity  Comparison  Between  Experiment,  KE,  and  MV  Sims  After 
Perforating  5083  Aluminum  Targets. 

residual  velocities  for  the  KE  Sims  are  closer  to  the  experimentally  determined  residual  velocities 
than  are  the  MV  Sims,  for  all  but  one  datum. 


3.3  BLINT  Model  Parameters.  To  show  the  differences  that  using  the  BLINT  model  can  make 
in  predictions,  three  of  the  KE  Sims  for  the  5083  aluminum  targets  were  repeated  without  the  BUNT 
model  active  (KE-NB-NC-P  Sims).  The  repeated  simulations  include  the  hemi-nose  penetrator 
shown  in  Figure  2b  (which  remained  rigid),  the  hemi-nose  penetrator  shown  in  Figure  2d  (which  was 
visibly  deformed/eroded),  and  the  ogival-nose  penetrator  at  the  highest  velocity  tested  (Figure  5e, 
which  remained  rigid).  Figure  1 1  shows  the  predicted  shapes  of  the  three  KE  Sims  that  were 
repeated  without  the  BUNT  model  active.  The  hemi-nose  penetrator  that  was  previously  predicted 
to  remain  rigid  (Figure  2b)  is  now  excessively  deformed  (Figure  11a).  The  hemi-nose  penetrator  that 
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Figure  10.  Residual  Velocity  Comparison  Between  Experiment,  KE,  and  MV  Sims  After 
Perforating  7039  Aluminum  Targets. 

was  previously  predicted  to  erode  (Figure  2b)  is  now  predicted  to  deform  and  erode  excessively 
(Figure  1  lb).  Finally,  with  the  BLINT  model  active,  all  the  ogival-nose  penetrators  were  predicted 
to  remain  rigid  at  all  striking  velocities  simulated;  however,  without  the  BLINT  model,  erosion  is 
predicted  (Figure  1  Ic).  Figure  12  compares  the  residual  velocities  of  the  KE  Sims  with  and  without 
the  BLINT  model  to  the  experimentally  determined  values,  hi  all  cases,  the  predicted  residual 
velocities  for  the  KE-NB-NC-P  Sims  underpredicted  those  with  the  BLINT  model  active.  The 
degree  to  which  the  KE-NB-NC-P  Sims  underpredicted  experiment  seems  to  increase  at  the  lower 
initial  striking  velocities. 

To  examine  the  effect  the  yield  strength  correction  factor  has  on  simulation  results,  the  hemi-  and 
ogival-nose  KE  Sims  against  the  5083  aluminum  target  were  repeated  without  the  strength  correction 
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Figure  11.  Residual  Shapes  of  KE-NB-NC-P  Sims  After  Perforating  5083  Aluminum  Targets. 


factor  (KE-BL-NC-P  Sims).  Residual  shapes  and  plastic  strains  of  the  hemi-  and  ogival-nose 
penetrators  are  shown  in  Figures  13  and  14,  respectively.  The  hemi-nose  KE-BL-NC-P  Sims  (Figure 
13)  appear  to  deform  in  a  manner  similar  to  the  hemi-nose  KE-NB-NC-P  Sims  (Figures  11a  and 
1  lb).  The  baseline  KE  Sims  for  the  ogival-nose  penetrators  shown  in  Figures  4b-4e  remained  rigid. 
Without  the  strength  correction  factor,  the  ogival-nose  penetrator  shows  visible  deformation  at  the 
very  tip  of  the  nose  at  a  striking  velocity  of  1,286  m/s  (Figure  14b)  that  progressively  increases  with 
striking  velocity  (Figures  14c-14e).  The  residual  velocity  for  the  KE-BL-NC-P  Sims  is  compared 
to  the  experimental  results  and  the  KE-BL-CO-P  baseline  simulations  in  Figure  15.  The  KE-BL- 
NC-P  Sims  underpredict  the  experimental  determined  residual  velocity  to  a  greater  extent  then  the 
baseline  KE  Sims.  At  striking  velocities  of  1,286  m/s  and  1,399  m/s,  the  predicted  residual 
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Figure  12.  Comparison  of  Residual  Velocity  for  KE-NB-NC-P  Sims  With  Experiment  and 
Baseline  KE  Sims  After  Perforating  5083  Aluminum  Targets. 

velocities  for  the  ogival-nose  penetrator  are  close  to  the  baseline  predictions,  because  the  penetrator 
was  only  slightly  deformed  (Figures  14b  and  14c). 

3.4  Failure  Model  Effects.  In  examining  the  effects  of  nose  shape  on  the  threshold  velocity  at 
which  tungsten  alloy  penetrators  transition  from  rigid  body  to  eroding  rods  when  penetrating  deep 
aluminum  targets,  it  has  been  shown  that  the  failure  model  used  can  influence  the  predicted 
penetration  depth  [8,  9].  Magness  has  suggested  a  more  appropriate  failure  model  would  be  one 
based  on  strain  [26],  the  reasons  for  which  are  given  in  Magness  [27].  Therefore,  all  KE  Sims  for 
the  5083  aluminum  targets  using  the  BLINT  model  were  redone  using  the  Johnson-Cook  failure  with 
all  but  the  first  parameter  set  to  zero,  such  that  all  materials  would  fail  when  they  exceeded  a 
threshold  value  of  150%  strain  or  when  they  exceeded  the  threshold  value  of  tensile  pressure 
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Figure  13.  Residual  Shapes  for  Hemi-Nose  Penetrator  KE-BL-NC-P  Sims  After  Perforating 
5083  Aluminum  Targets. 


reported  earlier.  The  predicted  shapes  for  the  KE-BL-CO-S  hemi-nose  penetrators  (Figure  16)  do 
not  differ  significantly  from  the  baseline  KiE  Sims  (Figures  2b-2d),  and  the  predicted  threshold 
velocity  remains  the  same  (between  1,147  and  1,201  m/s).  Even  the  levels  of  effective  plastic  strain 
observed  are  similar.  The  KE-BL-CO-S  ogival-nose  penetrators  (Figure  17)  remained  rigid  for  all 
striking  velocities  considered  as  did  the  baseline  KE  Sim  counterparts  (Figures  4b-4e);  however  the 
KE-BL-CO-S  experienced  larger  plastic  strains  for  the  same  striking  velocity.  Using  a  strain-based 
failure  criteria  reduced  the  predicted  residual  velocity  for  all  simulations  when  compared  to  the 
baseline  KE  Sims  and  experimental  results  other  than  for  the  hemi-nose  penetrator  with  a  striking 
velocity  of  1,201  m/s  (Figure  18). 
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Figure  14.  Residual  Shapes  for  Ogival-Nose  Penetrator  KE-BL-NC-P  Sims  After  Perforating 
5083  Aluminum  Targets. 


4.  Conclusions 


It  is  known  that  the  constitutive  response  of  tungsten  alloy  is  dependent  on  strain,  strain  rate, 
temperature,  percentage  tungsten  content,  tungsten  grain  size,  and  amount  of  swaging  [28].  In 
addition,  for  solid-solid  impacts  at  velocities  of  500-2,000  m/s,  impact  pressures  rapidly  decay  to 
values  comparable  to  the  strength  of  the  material;  therefore,  the  constitutive  model  is  of  primary 


striking  Velocity  (m/s) 

Figure  15.  Comparison  of  Residual  Velocity  for  KE-BL-NC-P  Sims  With  Experiment  and 
Baseline  KE  Sims  After  Perforating  5083  Aluminum  Targets. 

importance  and  the  EOS  is  of  secondary  importance  [29].  Therefore,  it  is  unrealistic  to  expect  the 
results  of  the  simulations,  with  the  constitutive  model  approximations  used  for  the  95W-2.5Ni- 
1.0Fe-L5Co,  21%  swaged  alloy,  to  provide  an  exact  match  with  the  experimental  data. 
Nevertheless,  the  following  conclusions  are  offered. 

The  BLINT  model  represents  a  drastic  improvement  in  the  predictive  capabilities  in  the  CTH 
hydrocode  for  certain  types  of  penetration  scenarios,  such  as  rigid  body  penetrations.  The  hydrocode 
was  able  to  predict  the  effect  of  rod  nose-shape  on  the  threshold  velocity  at  which  transition  from 
rigid  body  to  eroding  rod  penetration  occurs.  The  code  successfully  estimated  the  velocity  at  which 
the  hemi-nose  penetrators  begin  to  deform  in  penetrating  either  5083  or  7039  aluminum.  The 
transition  occurs  at  a  much  higher  striking  velocity  for  an  ogival-nose  penetrator  than  for  a  hemi- 
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Figure  16.  Residual  Shapes  for  Hemi-Nose  Penetrator  KE-BL-CO-S  Sims  After  Perforating 
5083  Aluminum  Targets. 


nose  penetrator.  Experimentally,  the  transition  occurred  with  the  onset  of  plastic  deformation  of  the 
hemi-nose  penetrator  at  a  striking  velocity  of  1,201  m/s  to  full  erosion  at  a  striking  velocity  of 
1,283  m/s  when  impacting  a  5083  aluminum  target.  The  KE  Sims  predicted  the  transition  velocity 
for  the  hemi-nose  penetrator  would  occur  between  a  striking  velocity  of  1,147  and  1,201  m/s  for  the 
5083  aluminum  target.  The  MV  Sims  hemi-nose  penetrator  only  predicted  minimal  plastic 
deformation  in  the  range  of  striking  velocities  between  1,147  and  1,283  m/s  for  the  5083  aluminum 
target.  Experimentally,  the  transition  velocity  for  the  ogival-nose  penetrator  impacting  the  5083 
aluminum  target  was  never  determined.  The  KE  Sims  predicted  the  visible  plastic  deformation 
would  occur  at  a  striking  velocity  between  1,900  and  2,000  m/s  and  erosion  of  the  penetrator  would 
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Figure  17.  Residual  Shapes  for  Ogival-Nose  Penetrator  KE-BL-CO-S  Sims  After  Perforating 
5083  Aluminum  Targets. 


occur  at  a  striking  velocity  between  2,000  and  2,100  m/s  for  the  ogival-nose  penetrator  impacting 
a  5083  aluminum  target.  The  transition  velocity  was  not  determined  for  the  ogival-nose  MV  Sims. 

Experimentally,  for  the  hemi-nose  penetrator  impacting  a  7039  aluminum  target,  plastic 
deformation  occurred  between  a  striking  velocity  of  1,038  and  1,093  m/s  and  erosion  occurred 
between  a  striking  velocity  of  1,093  and  1,198  m/s.  The  same  result  was  predicted  by  the  hemi-nose 


Figure  18.  Comparison  of  Residual  Velocity  for  KE-BL-CO-S  Sims  With  Experiment  and 
Baseline  KE  Sims  After  Perforating  5083  Aluminum  Targets. 

KE  Sims  for  the  7039  aluminum  target.  The  MV  Sims  predicted  visible  plastic  deformation  would 
occur  between  a  striking  velocity  of  1,093  and  1,198  m/s  for  a  7039  aluminum  target. 
Experimentally,  the  transition  velocity  occurred  between  a  striking  velocity  of  1,755  and  1,768  m/s 
for  the  ogival-nose  penetrator  impacting  a  7039  aluminum  target.  The  ogival-nose  KE  Sims 
predicted  visible  plastic  deformation  would  occur  at  a  striking  velocity  between  1,800  and  1,900  m/s 
and  erosion  at  a  striking  velocity  between  1,900  and  2,000  m/s  when  impacting  a  7039  aluminum 
target. 

While  the  simulations  were  not  in  complete  agreement  with  experiment,  they  did  a  reasonable 
job  at  predicting  the  transition  velocities  and  did  an  excellent  job  predicting  all  the  experimental 
trends.  The  simulations  predicted  the  large  difference  in  the  transition  velocity  for  a  hemi-nose  and 
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ogival-nose  penetrator,  and  they  correctly  predicted  the  effect  of  target  materials.  In  general,  the 
baseline  KE  Sims  did  much  better  at  predicting  final  rod  shape  than  did  the  MV  Sims. 

The  baseline  KE  Sims  did  better  in  predicting  residual  velocity  than  did  the  MV  Sims  for  both 
target  materials.  Both  the  ogival-nose  KE  and  MV  Sims  were  witiiin  8.1%  of  the  experimentally 
determined  residual  velocity,  except  for  one  datum.  The  hemi-nose  MV  Sims  were  closer  to 
predicting  the  experimentally  determined  residual  velocity  for  die  5083  aluminum  target,  and  the  KE 
Sims  were  closer  to  predicting  the  experimentally  determined  residual  velocity  for  the  7039 
aluminum  target.  For  the  5083  aluminum  target,  the  hemi-nose  KE  and  MV  Sims  predicted  the 
residual  velocity  within  9%.  For  the  7039  aluminum  target,  the  hemi-nose  KE  and  MV  Sims  under 
predicted  residual  velocity  by  as  much  as  24%. 

It  has  been  shown  that  without  the  BUNT  model  CTH  cannot  accurately  model  transition 
velocity  or  the  perforation  of  finite  target  plates.  Hemi-  and  ogival-nose  penetrators  that  were 
previously  predicted  to  remain  rigid,  eroded  and  deformed  excessively  without  the  BUNT  model 
active.  A  hemi-nose  penetrator  previously  predicted  to  erode,  eroded  and  deformed  excessively. 
Due  to  excessive  deformation  predicted  with  the  BUNT  model  inactive,  the  predicted  residual 
velocities  were  significantly  underpredicted.  It  was  also  shown  that  the  yield  strength  correction 
factor  greatly  influenced  simulation  results.  Without  the  yield  strength  correction  factor,  hemi-nose 
penetrator  simulation  results  were  similar  to  those  without  the  BUNT  model  active.  Without  the 
yield  strength  correction  factor,  the  ogival-nose  penetrators  deformed  and  eroded  prematurely.  As 
a  result,  most  simulations  run  without  the  yield  strength  correction  factor  under  predicted  residual 
velocity  by  a  significant  amount. 

Finally,  it  was  shown  that  the  choice  of  material  failure  models  could  influence  simulation 
results.  The  residual  shapes  and  residual  velocities  predicted  using  a  strain-based  failure  model  did 
not  differ  significantly  for  the  hemi-nose  simulations.  The  predicted  residual  velocity  was  greatly 
influenced  by  the  choice  of  failure  models.  No  attempt  was  made  to  predict  the  transition  velocity 
for  ogival-nose  simulations  using  a  strain-based  failure  model.  However,  the  effective  plastic  strains 
in  the  ogival-nose  penetrators  were  greater  than  the  baseline  simulations,  suggesting  that  the 


27 


transition  velocity  might  be  lower.  The  residual  velocities  predicted  by  the  ogival-nose  penetrator 
using  a  strain-based  failure  model  were  somewhat  reduced. 
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Appendix  A: 

Input  for  Hemi-Nose  KE-BL-CO-P  Sims  vs.  5083  Aluminum 
Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes 
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*  Run  History 

* 

*eor*  cgenin 
* 

nose  shape  tests:  hemi  v-1 147  target-5083  al  finite 
* 

control 

ep 

mmp 

viscosity  bl-.l  bq=2  bs»0.1 

endcontrol 

* 

mesh 

block  1  geom-2dc  type-e 
xO  0.00 

xl  n  40  w  1.6891  rat  1.0 
x2  n  45  dxf  0.0422275  rat  1 .05 
endx 

yO  -10.219055 

yl  n  832  dyf  0.0422275  rat  1.0 
endy 

xactive  0.0  0.4 

yactive  -10.219055  0.0 
endblock 

endmesh 

* 

insertion 
block  1 

package  ‘5083  Al  target’ 
material  1 
numsub  50 
insert  box 
pl=0.0  0.0 

p2-7.6  7.62 

endinsert 
endpackage 

package  ‘w  alloy  rod  nose’ 
material  2 
numsub  100 

* 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

* 

yvel  1.147e5 
insert  circle 
ce-0.0  -0.33782 


35 


r- 0.33782 
endinsert 
endpackage 

package  ‘w  alloy  rod  body’ 
material  2 
numsub  50 

* 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

* 

yvel  1.147e5 
insert  box 
pl-0.0  -0.33782 

p2-0.33782 -10.1346 
endinsert 
endpackage 
endblock 

endinsertion 

* 

epdata 

vpsave 

♦ 

*  NOTE:  poisson  ratio  from  for  Tungsten  from  Metals  Handbook,  that  for  5083 

*  Aluminum  from  Forrestal,  M.J.,  V.  K.  Luk,  and  N.  S.  Brar,  “Perforation  of 

*  aluminum  with  conical-nose  projectiles’’,  Mechanic  and  Materials  10  (1990) 

*  pp.  97-105. 

* 

*  NOTE:  properties  for  5083  aluminum  from  Silling,  S.  A.,  “(TTH  Reference  Manual: 

*  Boundary  Layer  Algorithm  for  Sliding  Interfaces  in  Two  Dimensions”,  Sandia 

*  National  Laboratories  Report  SAND93-2487,  January  1994. 

* 

*  NOTE:  5083  aluminum  using  undocumented  power  law 

* 

matep  1 

johnson-cook«‘USER  ’ 

ajo=-2.76e9  bJo=254.7  cjo=0.0 
mjo=1.0  njo-0.084  tjo«=6.68e-2 

poisson  0.333 

* 

*  NOTE:  actual  tungsten  alloy  was  95W-2.5Ni-l.0Fe-l.5Co  (21%  swaged)  with  rO-18.1. 

*  95W-3.5Ni- 1 .5Fe  is  being  used  to  approximate  the  w  alloy. 

* 

matep  2 

steinberg-‘TUNGSTEN_NI_FE’ 

r0st=18.16  tm0st=0. 195002  atmst=1.3 
gmOst-1.67  ast-1.03e-12  bst-1.76396 
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nst-0.13  clst-0.0  0281-0.0 
g0st-1.45el2  btst-7.7  eist-0.0 
ypst-0.0  ukst-0.0  ysmst-0.0 
yast-0.0  y0st-18.7e9  ymst=40.e9 
poisson-0.280 

* 

*  NOTE:  parameters  for  boundary  layer  algorithm  taken  similar  to: 

*  Kmetyk,  L.  N.  and  P.  Yarrington,  “CTH  Analysis  of  Steel  Rod  Penetration 

*  Into  Aluminum  and  Concrete  Targets  with  Comparisons  to  Experimental  Data”, 

*  Sandia  National  Laboratories  Report  SAND94-1498,  October  1994, 

* 

*  NOTE:  if  no  BLINT  model,  next  line  is  commented  out  If  no  yield  strength  correction 

*  factor  ‘corr’  is  omitted 

* 

blint  1  soft  1  hard  2  wsl  0,084455  wbl  0.084455  fnc-0.0  corr 
mix  3 

endq) 

* 

tracer 
block  1 

addO.  0.  to  0. -10.134  n  10 
endtracer 
♦ 

edit 
block  1 
noexpanded 
endblock 

endedit 

♦ 

eos 

* 

*  NOTE:  EOS  properties  from  cth  mgrun  library 

* 

*  NOTE:  5083-H131  Aluminum  eos  approximated  with  6061-t6  Aluminum 

*  density  reduced  to  reflect  that  for  5083  Aluminum  from 

*  Metal’s  Handbook. 

* 

matl  mgrun  eos-6061-t6_al  1O-2.66  cs=0.534e6  s-1.4 
gO-1.97  cv-1.07ell 

mat2mgruneos-tungsten_nir0-18.16  cs-0.403e6  s-1.237 
gO-1.67  cv-1.66el0 

endeos 

* 

*eor*  cthin 
♦ 
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nose  shape  tests:  hemi  v-1 147  target-5083  al  finite 
* 

control 
tstop  300.e-6 
rdumpf  3600 
cpshift  999. 

endcontrol 

* 

restart 

*  file-‘rsctr 
time-O.e-6 

endr 

♦ 

cellthermo 

mmp 

ntbad-99999 

endc 

* 

convct 

* 

*  NOTE:  if  KE  Sim  convection— 0,  if  MV  Sim  convection— 1 

* 

convection-0 

interface=high_resolution 

endconvct 

* 

edit 

shortt 

time  0.  dtffequency  150.e-6 
ends 
longt 

time  0.  dtffequency  600.e-6 
endl 
plott 

time  0.  dtfrequency  50.e-6 
endp 
histt 

time  0.  dtffequency  0.3e-6 

htracerl 

htracer2 

htracer3 

htracer4 

htracerS 

htracerd 

htracer? 
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htracerS 

htracer9 

htracerlO 

endh 

ende 

♦ 

boundary 
bhydro 
block  1 
bxbot  0 
bxtop  2 
bybot2 
bytop  2 
endb 
endh 

endb 

♦ 

fracts 
pressure 
pfracl  -4.5e9 
pfrac2  -35.0e9 
pfinix  -1.0e20 
pfvoid  -1.0e20 

endf 

♦ 

*eor*  pltinp 


Intentionally  left  blank. 
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Appendix  B: 


Input  for  Ogival-Nose  KE-BL-CO-P  Sims  vs.  5083  Aluminum 
Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes 


Intentionally  1£ft  blank. 
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**  Run  History 
* 

*eor*  cgenin 
* 

nose  shape  tests:  ogive  v-1286  target-5083  al  finite 
* 

control 

ep 

nunp 

viscosity  bl-.l  bq-2  bs-0.1 

endcontrol 

* 

mesh 

block  1  geom=2dc  type=e 
xOO.OO 

xl  n  40  w  1.6891  rat  1.0 
x2n45  dxf  0.0422275  rat  1.05 
endx 

yO -10.219055 

yl  n  832  dyf  0.0422275  rat  1.0 
endy 

xactive  0.0  0.4 

yactive  -10.219055  0.0 
endblock 

endmesh 

♦ 

insertion 
block  1 

package  ‘5083  Al  target’ 
material  1 
numsub  50 
insert  box 
pl-0.0  0.0 

p2-7.6  7.62 

endinsert 
endpackage 

package  ‘w  alloy  ogive  nose  rod’ 
material  2 
numsub  100 

♦ 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 


yvel 1.286e5 
insert  uds 
point  0.0  0.0 
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point  0.050437038  -0.08003015622 
point  0.095867718  -0.1600603124 
point  0.136665799  -0.2400904687 
point  0.173135779  -0.3201206249 
point  0.205527612  -0.400150781 1 
point  0.234047376  -0.4801809373 
point  0.258865148  -0.56021 10935 
point  0.280120900  -0.6402412498 
point  0.297928948  -0.7202714060 
point  0.312381319  -0.8003015622 
point  0.323550306  -0.8803317184 
point  0.331490355  -0.9603618746 
point  0.336239443 -1.040392031 
point  0.337820  -1.120422187 
point  0.337820  -10.1346 
point  0.0  -10.1346 

point  0.0  0.0 

endinsert 
endpackage 
endblock 

endinsertion 

* 

epdata 

vpsave 

* 

*  NOTE:  poisson  ratio  from  for  Tungsten  from  Metals  Handbook,  that  for  5083 

*  Aluminum  from  Foirestal,  M.J.,  V.  K.  Luk,  and  N.  S.  Brar,  “Perforation  of 

*  aluminum  with  conical-nose  projectiles”.  Mechanic  and  Materials  10  (1990) 

*  pp.  97-105. 

* 

*  NOTE:  properties  for  5083  aluminum  from  Silling,  S.  A.,  “CTH  Reference  Manual: 

*  Boundary  Layer  Algorithm  for  Sliding  Interfaces  in  Two  Dimensions”,  Sandia 

*  National  Laboratories  Report  SAND93-2487,  January  1994. 

*  NOTE:  5083  aluminum  using  undocumented  power  law 

* 

matep  1 

johnson-cook^’USER’ 

ajo— 2.76e9  bjo««254.7  cjo-0.0 
mjo-1.0  njo-0.084  tjo-6.68e-2 

poisson  0.333 

* 

*  NOTE:  actual  tungsten  alloy  was  95W-2.5Ni-l.0Fe-l.5Co  (21%  swaged)  with  i0=18.1. 

*  95W-3.5Ni- 1 .5Fe  is  being  used  to  approximate  the  w  alloy. 
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matep2 

steinberg-‘TUNGSTEN_NI_FE’ 

rOst-18.16  tmOst-0. 195002  atmst-1.3 
gmOst-1.67  ast“1.03e-12  bst>=1.76396 
nst-0.13  clst-0.0  c2st«0.0 

gOst- 1 .45e  1 2  btst-7.7  eist-0.0 

ypst-0.0  ukst-0.0  ysmst-0.0 

yast-=0.0  y0st=18.7e9  ymst=40.e9 

poisson-0.280 

* 

*  NOTE:  parameters  for  boundary  layer  algorithm  taken  similar  to: 

*  Kmetyk,  L.  N.  and  P.  Yarrington,  “CTH  Analysis  of  Steel  Rod  Penetration 

*  Into  Aluminum  and  Concrete  Targets  with  Comparisons  to  Experimental  Data”, 

*  Sandia  National  Laboratories  Report  SAND94-1498,  October  1994. 

* 

*  NOTE:  if  no  BUNT  model,  next  line  is  commented  out.  If  no  yield  strength  correction 

*  factor  ‘corr’  is  omitted 

* 

blint  1  soft  1  hard  2  wsl  0.084455  wbl  0.084455  fric-O.O  corr 
mix  3 

endep 

* 

tracer 
block  1 

add  0.0.  to  0.  -10.1345  n  10 

endtracer 

* 

edit 
block  1 
noexpanded 
endblock 

endedit 

♦ 

eos 

* 

*  NOTE:  EOS  properties  from  cth  mgmn  library 

♦ 

*  NOTE:  5083-H131  Aluminum  eos  approximated  with  6061-t6  Aluminum 

*  density  reduced  to  reflect  that  for  5083  Aluminum  from 

*  Metal’s  Handbook. 

* 

matl  mgrun  eos“6061-t6_al  r0=2.66  cs-0.534e6  s=1.4 
g0=1.97  cv-1.07ell 

mat2mgruneos“tungsten_nir0-18.16  cs-0.403e6  s=1.237 
g0«1.67  cv-1.66el0 
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endeos 

* 

*eor*  cthin 
* 

nose  shape  tests:  ogive  v-1286  target-5083  al  finite 
* 

control 
tstop300.e-6 
rdumpf  3600 
cpshift  999. 

endcontrol 

* 

restart 

*  file-‘rsctr 
time-O.e-6 

endr 

* 

cellthermo 

mmp 

ntbad-99999 

endc 

* 

convct 

* 

*  NOTE:  if  KE  Sim  convection-0,  if  MV  Sim  convection- 1 

* 

convection-0 

interface-high_resolution 

endconvct 

* 

edit 

shortt 

time  0.  dtffequency  150.e-6 
ends 
longt 

time  0.  dtfrequency  600.e-6 
endl 
plott 

time  0.  dtfi'equency  50.e-6 
endp 
histt 

time  0.  dtfi'equency  0.3e-6 

htracerl 

htracer2 

htracer3 
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htracer4 

htracerS 

htracer6 

htracer? 

htracerS 

htracer9 

htracerlO 

endh 

ende 

* 

boundary 
bhydro 
block  1 
bxbot  0 
bxtop  2 
bybot  2 
bytop  2 
endb 
endh 

endb 

* 

fracts 
pressure 
pfracl  -4.5e9 
pfrac2  -35.0e9 
pfinix  -1.0e20 
pfvoid  -1.0e20 

endf 

* 

*eor*  pltinp 


Intentionally  left  blank. 
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Appendix  C: 

Input  for  Hemi-Nose  KE-BL-CO-P  Sims  vs.  7039  Aluminum 
Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes 
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Intentionally  left  blank. 


*  Run  History 

* 

*eor*  cgenin 
♦ 

nose  shape  tests:  hemi  v=1038  target-7039  al  finite 

control 

ep 

nrnip 

viscosity  bl-.l  bq-2  bs-0.1 

endcontrol 

* 

mesh 

block  1  geom-2dc  type-e 
xOO.OO 

xl  n  40  w  1.6891  rat  1.0 
x2n45  dxf  0.0422275  rat  1.05 
endx 

yO  -10.219055 

yl  n  832  dyf  0.0422275  rat  1.0 
endy 

xactive  0.0  0.4 

yactive  -10.219055  0.0 
endblock 

endmesh 

* 

insertion 
block  1 

package  ‘7039  Al  target’ 
material  1 
numsub  50 
insert  box 
pl-0.0  0.0 

p2=7.6  7.62 

endinsert 
endpackage 

package  ‘w  alloy  rod  nose’ 
material  2 
numsub  100 

♦ 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

* 

yvel  1.038e5 
insert  circle 
ce-0.0  -0.33782 
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r-  0.33782 
endinsert 
endpackage 

package  *w  alloy  rod  body’ 
material  2 
numsub  50 

* 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

* 

yvel 1.038e5 
insert  box 
pl-0.0  -0.33782 

P2-0.33782 -10.1346 
endinsert 
endpackage 
endblock 

endinsertion 

* 

epdata 

vpsave 

♦ 

*  NOTE:  poisson’s  ratio  from  Metals  Handbook 

* 

matep  1 

johnson-cook-‘7039_ALUMINUM’ 

ajo=3.3672e9  bjo=3.4293e9  cjo=0.01 
mjo-1.0  njo=.41  tjo-7.76342e-2 

poisson  0.345 

♦ 

*  NOTE:  actual  tungsten  alloy  was  95W-2.5Ni-l.0Fe-l.5Co  (21%  swaged)  with  rO-18.1. 

*  95 W-3 .5Ni- 1 .5Fe  is  being  used  to  approximate  the  w  alloy. 

* 

matep  2 

steinberg=‘TUNGSTEN_NI_FE’ 

r0st“18.16  tm0st“0. 195002  atmst— 1.3 
gm0st=1.67  ast=1.03e-12  bst-1.76396 
nst—0.13  clst—0.0  c2st— 0.0 

g0st-1.45el2  btst-7.7  eist-0.0 

ypst-0.0  ukst-0.0  ysmst-0.0 

yast-0.0  y0st“18.7e9  ymst-40.e9 
poisson=0.280 

* 

*  NOTE:  parameters  for  boundary  layer  algorithm  taken  similar  to: 

*  Kmetyk,  L.  N.  and  P.  Yarrington,  “CTH  Analysis  of  Steel  Rod  Penetration 

*  Into  Aluminum  and  Concrete  Targets  with  Comparisons  to  Experimental  Data”, 
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*  Sandia  National  Laboratories  Report  SAND94-1498,  October  1994. 

* 

*  NOTE:  if  no  BLINT  model,  next  line  is  commented  out.  If  no  yield  strength  correction 

*  factor  ‘corr’  is  omitted 

* 

blint  1  soft  1  hard  2  wsl  0.084455  wbl  0.084455  fiic-0.0  corr 
mix  3 

endep 

♦ 

tracer 
block  1 

add  0. 0.  to  0.  -10.134  n  10 

endtracer 

* 

edit 
block  1 
noexpanded 
endblock 

endedit 

* 

eos 

* 

*  NOTE:  EOS  properties  from  cth  mgrun  library 

* 

*  NOTE:  7039  Aluminum  eos  approximated  with  7075-t6  Aluminum 

*  density  reduced  to  reflect  that  for  7039  Aluminum  from 

*  Johnson  and  Cook  (1983). 

* 

matl  mgran  eos-7075-t6_al  r0=2.77  cs-0.520e6  s*=1.36 
gO-2.20  cv-1.07ell 

mat2  mgrun  eos=tungsten_ni  r0=  1 8. 1 6  cs=0.403e6  s=l  .237 
g0=1.67  cv-1.66el0 

endeos 

* 

*eor*  cthin 
♦ 

nose  shape  tests:  hemi  v-1038  target-7039  al  finite 
* 

control 
tstop  300.e-6 
rdumpf  3600 
cpshift  999. 

endcontrol 

* 

restart 


53 


*  file-‘rsctr 
time-O.e-6 

endr 

* 

cellthermo 

nunp 

ntbad-99999 

endc 

* 

convct 

* 

*  NOTE:  if  KE  Sim  convection-0,  if  MV  Sim  convection- 1 

* 

convection-0 

interface-high_resolution 

endconvct 

* 

edit 

shortt 

time  0.  dtffequency  150.e-6 
ends 
longt 

time  0.  dtffequency  600.e-6 
endl 
plott 

time  0.  dtffequency  50.e-6 
endp 
histt 

time  0.  dtfrequency  0.3e-6 
htracerl 
htracer2 
htracerS 
htracer4 
htracerS 
htracerti 
htracer? 
htracerS 
htracer9 
htracer  10 
endh 

ende 

♦ 

boundary 
bhydro 
block  1 
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bxbot  0 
bxtop  2 
bybot  2 
bytop  2 
endb 
endh 

endb 

* 

tracts 
pressure 
pfracl  -5.0e9 
pfrac2  -35.0e9 
pfmix  -1.0e20 
pfvoid  -1.0e20 

endf 

* 

*eor*  pltinp 


Intentionally  left  blank. 
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Appendix  D: 

Input  for  Ogival-Nose  KE-BL-CO-P  Sims  vs.  7039  Aluminum 
Targets  -  Changes  for  MV,  NB,  and  NC  Sims  Given  in  Notes 
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Intentionally  left  blank. 
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*  Run  History 

* 

*eor*  cgenin 
* 

nose  shape  tests:  ogive  v=1474  target=7039  al  finite 
* 

control 

ep 

inmp 

viscosity  bl=,l  bq=2  bs-=0.1 

endcontrol 

* 

mesh 

block  1  geom=2dc  type=e 
xO  0.00 

xl  n  40  w  1.6891  rat  1.0 
x2n45  dxf  0.0422275  rat  1.05 
endx 

yO  -10.219055 

yl  n  832  dyf  0.0422275  rat  1.0 
endy 

xactive  0.0  0.4 

yactive  -10.219055  0.0 
endblock 

endmesh 

* 

insertion 
block  1 

package  ‘7039  Al  target’ 
material  1 
numsub  50 
insert  box 
pl=0.0  0.0 

p2-7.6  7.62 

endinsert 
endpackage 

package  ‘ogive  nose  rod’ 
material  2 
numsub  100 

* 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

* 

yvel 1.474e5 
insert  uds 
point  0.0  0.0 
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point  0.050437038  -0.08003015622 
point  0.095867718  -0.1600603124 
point  0.136665799  -0.2400904687 
point  0.173135779  -0.3201206249 
point  0.205527612  -0.400150781 1 
point  0.234047376  -0.4801809373 
point  0.258865148  -0.56021 10935 
point  0.280120900  -0.6402412498 
point  0.297928948  -0.7202714060 
point  0.312381319  -0.8003015622 
point  0.323550306  -0.8803317184 
point  0.33 1490355  -0.960361 8746 
point  0.336239443  -1.040392031 
point  0.337820  -1 . 1204221 87 
point  0.337820  -10.1346 
point  0.0  -10.1346 

point  0.0  0.0 

endinsert 
endpackage 
endblock 

endinsertion 

* 

epdata 

vpsave 

* 

*  NOTE:  poisson’s  ratio  from  Metals  Handbook 

* 

matep  1 

johnson-cook-‘7039_ALUMB4UM’ 

ajo=3.3672e9  bjo=3.4293e9  cjo=0.01 
mjo-1.0  njo=.41  tjo=7.76342e-2 

poisson  0.345 

* 

*  NOTE:  actual  tungsten  aUoy  was  95W-2.5Ni-l.0Fe-l.5Co  (21%  swaged)  with  rO-18.1. 

*  95W-3.5Ni- 1 .5Fe  is  being  used  to  approximate  the  w  alloy. 

matep  2 

Steinberg-  ‘TUNGSTEN_NI_FE’ 

rOst-18.16  tmOst-0.195002  atmst-1.3 
gmOst-1.67  ast-1.03e-12  bst-1.76396 
nst-0.13  clst-0.0  c2st-0.0 
g0st-1.45el2  btst— 7.7  eist-0.0 
ypst-0.0  ukst-0.0  ysmst-0.0 
yast-0.0  y0st*=18.7e9  ymst-40.e9 
poisson-0.280 
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* 

*  NOTE:  parameters  for  boundary  layer  algorithm  taken  similar  to: 

*  Kmetyk,  L.  N.  and  P.  Yarrington,  “CTH  Analysis  of  Steel  Rod  Penetration 

*  Into  Aluminum  and  Concrete  Targets  with  Comparisons  to  Experimental  Data”, 

*  Sandia  National  Laboratories  Report  SAND94-1498,  October  1994. 

* 

*  NOTE:  if  no  BUNT  model,  next  line  is  commented  out.  If  no  yield  strength  correction 

*  factor  ‘corr’  is  omitted 

* 

blint  1  soft  1  hard  2  wsl  0.084455  wbl  0.084455  ftioO.O  corr 
mix  3 

endep 

* 

tracer 
block  1 

addO.  O.toO.  -10.1345  n  10 

endtracer 

* 

edit 
block  1 
noexpanded 
endblock 

endedit 

* 

eos 

♦ 

*  NOTE:  EOS  properties  from  cth  mgnm  library 

He 

*  NOTE:  7039  Aluminum  eos  approximated  with  7075-t6  Aluminum 

*  density  reduced  to  reflect  that  for  7039  Aluminum  from 

*  Johnson  and  Cook  (1983). 

He 

matl  mgnm  eos-7075-t6_al  i0-=2.77  cs-0.520e6  s-1.36 
g0=2.20  cv«1.07ell 

mat2mgruneos=tungsten_nir0-18.16  cs=0.403e6  s=1.237 
gO-1.67  cv»=1.66el0 

endeos 

He 

*eor*  cthin 

He 

nose  shape  tests:  ogive  v-1474  target-e7039  al  finite 

He 

control 
tstop  300.e-6 
rdumpf  3600 
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cpshift  999. 

endcontrol 

* 

restart 

*  file-«‘rsctr 
time-O.e-6 

endr 

* 

cellthermo 

mmp 

ntbad-99999 

endc 

iie 

convct 

* 

*  NOTE;  if  KE  Sim  convection-0,  if  MV  Sim  convection- 1 

* 

convection-0 

interface-high_resolution 

endconvct 

* 

edit 

shortt 

time  0.  dtfrequency  150.e-6 
ends 
longt 

time  0.  dtfrequency  600.e-6 
endl 
plott 

time  0.  dtfrequency  50.e-6 
endp 
histt 

time  0.  dtfrequency  0.3e-6 
htracerl 
htracer2 
htracerS 
htracer4 
htracerS 
htracerS 
htracer? 
htracerS 
htracer9 
htracerl  0 
endh 
ende 
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boundary 
bhydro 
block  1 
bxbot  0 
bxtop  2 
bybot  2 
bytop  2 
endb 
endh 

endb 

* 

fracts 
pressure 
pfracl  -5.0e9 
pfrac2  -35.0e9 
pfmix  -L0e20 
pfvoid  -1.0e20 
endf 


*eor*  pltinp 


Intentionally  left  blank. 
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Appendix  £: 

Input  for  Hemi-Nose  KE-BL-CO-S 
Sims  vs.  5083  Aluminum  Targets 


65 


Intentionally  left  blank. 


66 


*  Run  Histoiy 

* 

♦eor*  cgenin 
* 

nose  shape  tests:  hemi  v=l  147  target=5083  al  finite 

control 

ep 

mmp 

viscosity  bl=.l  bq-»2  bs-0.1 

endcontrol 

♦ 

mesh 

block  1  geom-=2dc  type=e 
xOO.OO 

xl  n  40  w  1.6891  rat  1.0 
x2n45  dxf 0.0422275  rat  1.05 
endx 

yO  -10.219055 

yl  n  832  dyf  0.0422275  rat  1.0 
endy 

xactive  0.0  0.4 

yactive  -10.219055  0.0 
endblock 

endmesh 

* 

insertion 
block  1 

package  ‘5083  Al  target’ 
material  1 
numsub  50 
insert  box 
pl-0.0  0.0 

p2=7.6  7.62 

endinsert 
endpackage 

package  ‘w  alloy  rod  nose’ 
material  2 
numsub  100 

* 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

♦ 

yvel  1.147e5 
insert  circle 
ce-0.0  -0.33782 
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r- 0.33782 
endinsert 
endpackage 

package  ‘w  alloy  rod  body’ 
material  2 
numsub  50 

* 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

* 

yvel  1.147e5 
insert  box 
pl-0.0  -0.33782 

p2-0.33782  -10.1346 
endinsert 
endpackage 
endblock 

endinsertion 

* 

epdata 

vpsave 


*  NOTE:  poisson  ratio  from  for  Tungsten  from  Metals  Handbook,  that  for  5083 

*  Aluminum  from  Foirestal,  M.J.,  V.  K.  Luk,  and  N.  S.  Biar,  “Perforation  of 

*  aluminum  with  conical-nose  projectiles”.  Mechanic  and  Materials  10  (1990) 

*  pp.  97-105. 


*  NOTE:  properties  for  5083  aluminum  from  Silling,  S.  A.,  “CTH  Reference  Manual: 

Boundary  Layer  Algorithm  for  Sliding  Literfaces  in  Two  Dimensions”,  Sandia 
National  Laboratories  Report  SAND93-2487,  January  1994. 


♦ 

♦ 

♦ 

*  NOTE:  5083  aluminum  using  undocumented  power  law 

♦ 

matep  1 

johnson-cook=‘USER’ 

ajo—2.76e9  bjo-254.7  cjo-0.0 
mjo-1.0  njo-0.084  tjo-6.68e-2 

poisson  0.333 

* 

*  NOTE:  for  strain  failure  at  150%  all  parameter  for  Johnson-Cook  fracture  model  except 

*  the  first  are  set  to  zero. 

* 

jfrac-‘USER’  jfdl=1.5  jfd2-0.0  jfd3-0.0  jfd4-0.0  jfd5-0.0 
jftm-O.Oj^fO— 4.5e9 

♦ 

*  NOTE:  actual  tungsten  alloy  was  95W-2.5Ni-l.0Fe-l.5Co  (21%  swaged)  with  i0=18.1. 
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*  95W-3.5Ni- 1 .5Fe  is  being  used  to  approximate  the  w  alloy. 

♦ 

matep  2 

steinberg=‘TUNGSTEN_NI_FE’ 

rOst-18.16  tmOst-0. 195002  atmst-1.3 
gmOst-1.67  ast-1.03e-12  bst-1.76396 
nst-0,13  clst-0.0  c2st-0.0 

gOst- 1 .45e  1 2  btst=7.7  eist=0.0 

ypst-4),0  ukst-0.0  ysmst-0.0 

yast-0.0  y0st-18.7e9  ymst-40.e9 
poisson-0.280 

* 

*  NOTE:  for  strain  failure  at  150%  all  parameter  for  Johnson-Cook  fracture  model  except 

*  the  first  are  set  to  zero. 

* 

jffac-‘USER’  jfdl-1 .5  jfd2-0.0  jfd3-0.0  jfd4-0.0  jfd5-0.0 
jftm-0.0  jl^fO— 35.e9 

* 

*  NOTE:  parameters  for  boundary  layer  algorithm  taken  similar  to: 

*  Kmetyk,  L.  N.  and  P.  Yarrington,  “CTH  Analysis  of  Steel  Rod  Penetration 

*  Into  Aluminum  and  Concrete  Targets  with  Comparisons  to  Experimental  Data”, 

*  Sandia  National  Laboratories  Report  SAND94-1498,  October  1994. 

* 

*  NOTE:  if  no  BUNT  model,  next  line  is  commented  out.  If  no  yield  strength  correction 

*  factor  ‘coir’  is  omitted 

* 

blint  1  soft  1  hard  2  wsl  0.084455  wbl  0.084455  fnc=0.0  corr 
mix  3 

endep 

♦ 

tracer 
block  1 

addO.  0.  to  0. -10.134  n  10 

endtracer 

* 

edit 
block  1 
noexpanded 
endblock 

endedit 

* 

eos 

* 

*  NOTE:  EOS  properties  from  cth  mgrun  library 

* 
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*  NOTE:  5083-H131  Aluminum  eos  approximated  with  6061-t6  Aluminum 

*  density  reduced  to  reflect  that  for  5083  Aluminum  from 

*  Metal’s  Handbook. 

* 

matl  mgruneos-6061-t6_al  r0<=2.66  cs-0.534e6  s-1.4 
gO-1.97  cv-1.07ell 

mat2mgruneos-=tungsten_ni  10=18.16  cs=0.403e6  s=1.237 
gO-1.67  cv-1.66el0 

endeos 

* 

*eor*  cthin 
* 

nose  shape  tests:  hemi  v=1147  target=5083  al  finite 
♦ 

control 
tstop  300.e-6 
rdumpf  3600 
cpshift  999. 

endcontrol 

* 

restart 

*  file=‘rsctr 
time-O.e-6 

endr 

* 

cellthermo 

mmp 

ntbad=99999 

endc 

% 

convct 

* 

*  NOTE:  if  KE  Sim  convection=0,  if  MV  Sim  convection=l 

* 

convection-0 

interface-high_resolution 

endconvct 

* 

edit 

shortt 

time  0.  dtfrequency  150.e-6 
ends 
longt 

time  0.  dtfrequency  600.e-6 
endl 
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plott 

time  0.  dtfrequency  50.e-6 
endp 
histt 

time  0.  dtfrequency  0.3e-6 
htracerl 
htracer2 
htracerS 
htracer4 
htracerS 
htracer6 
htracer? 
htracerS 
htracer9 
htracerl  0 
endh 

ende 

* 

boundary 
bhydro 
block  1 
bxbot  0 
bxtop  2 
bybot  2 
bytop  2 
endb 
endh 

endb 

* 

fracts 
pressure 
pfracl  -4.5e9 
pfrac2  -35.0e9 
pfmix  -1.0e20 
pfvoid  -1.0e20 
endf 
* 

*eor*  pltinp 


Intentionally  left  blank. 
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Appendix  F: 

Input  for  Ogival-Nose  KE-BL-CO-S 
Sims  vs.  5083  Aluminum  Targets 


73 


Intentionally  left  blank. 
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*  Run  History 

♦ 

*eor*  cgenin 
♦ 

nose  shape  tests;  ogive  v-=1286  target=5083  al  finite 
* 

control 

ep 

mmp 

viscosity  bl-.l  bq-2  bs-0.1 

endcontrol 

* 

mesh 

block  1  geom=2dc  type=e 
xOO.OO 

xl  n  40  w  1.6891  rat  1.0 
x2n45  dxf  0.0422275  rat  1.05 
endx 

yO  -10.219055 

yl  n  832  dyf  0.0422275  rat  1.0 
endy 

xactive  0.0  0.4 

yactive  -10.219055  0.0 
endblock 

endmesh 

* 

insertion 
block  1 

package  ‘5083  Al  target’ 
material  1 
numsub  50 
insert  box 
pl=0.0  0.0 

p2=7.6  7.62 

endinsert 
endpackage 

package  ‘w  alloy  ogive  nose  rod’ 
material  2 
numsub  100 

* 

*  NOTE:  striking  velocity  (yvel)  is  changed  below 

* 

yvel  1.286e5 
insert  uds 
point  0.0  0.0 
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point  0.050437038  -0.08003015622 
point  0.095867718  -0.1600603124 
point  0.136665799  -0.2400904687 
point  0.173135779  -0.3201206249 
point  0.205527612  -0.400150781 1 
point  0.234047376  -0.4801809373 
point  0.258865148  -0.5602110935 
point  0.280120900  -0.6402412498 
point  0.297928948  -0.7202714060 
point  0.312381319  -0.8003015622 
point  0.323550306  -0.8803317184 
point  0.331490355  -0.9603618746 
point  0.336239443  -1.040392031 
point  0.337820  -1.120422187 
point  0.337820  -10.1346 
point  0.0  -10.1346 

point  0.0  0.0 

endinsert 
endpackage 
endblock 

endinsertion 

* 

epdata 

vpsave 

♦ 

*  NOTE:  poisson  ratio  from  for  Tungsten  from  Metals  Handbook,  that  for  5083 

*  Aluminum  from  Forrestal,  M.J.,  V.  K.  Luk,  and  N.  S.  Brar,  “Perforation  of 

*  aluminum  with  conical-nose  projectiles”.  Mechanic  and  Materials  10  (1990) 

*  pp.  97-105. 

* 

*  NOTE:  properties  for  5083  aluminum  from  Silling,  S.  A.,  “CTH  Reference  Manual: 

*  Boundary  Layer  Algorithm  for  Sliding  Interfaces  in  Two  Dimensions”,  Sandia 

*  National  Laboratories  Report  S AND93-2487,  January  1 994. 

* 

*  NOTE:  5083  aluminum  using  undocumented  power  law 

♦ 

matep  1 

johnson-cook-=‘USER’ 

ajo— 2.76e9  bjo— 254.7  cjo-0.0 
mjo-1.0  njo=0.084  tjo-=6.68e-2 

poisson  0.333 

* 

*  NOTE:  for  strain  failure  at  150%  all  parameter  for  Johnson-Cook  fracture  model  except 

*  the  first  are  set  to  zero. 
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jfrac-‘USER’  jfdl-1 .5  jfd2-0.0  jfd3-0.0  jfd4-0.0  jfd5-0.0 
jftm-0.0  jfpfO— 4.5e9 

* 

*  NOTE:  actual  tungsten  alloy  was  95W-2.5Ni-l.0Fe-l.5Co  (21%  swaged)  with  rO-18.1. 

*  95W-3.5Ni- 1 .5Fe  is  being  used  to  approximate  the  w  alloy. 

* 

matep  2 

steinberg-‘TUNGSTEN_NI_FE’ 

iOst-18.16  tmOst-0. 195002  atmst-1.3 
gm0st=1.67  ast-1.03e-12  bst-1.76396 
nst-0.13  clst-0.0  c2st-0.0 
gOst^l  .45el  2  btst=7.7  eist-0.0 
ypst-0.0  ukst-0.0  ysmst-0.0 
yast=0.0  y0st*=18.7e9  ymst=40.e9 
poisson-0.280 

* 

*  NOTE:  for  strain  failure  at  150%  all  parameter  for  Johnson-Cook  fracture  model  except 

*  the  first  are  set  to  zero. 

* 

jfrac-‘USER’  jfdl-1.5  jfd2-0.0  jfd3-0.0  jfd4-0.0  jfd5-0.0 
jftm-=0.0  jfpf0=-35.e9 

* 

*  NOTE:  parameters  for  boundary'  layer  algorithm  taken  similar  to: 

*  Kmetyk,  L.  N.  and  P.  Yanington,  “CTH  Analysis  of  Steel  Rod  Penetration 

*  Into  Aluminum  and  Concrete  Targets  with  Comparisons  to  Experimental  Data”, 

*  Sandia  National  Laboratories  Report  SAND94-1498,  October  1994. 

* 

*  NOTE:  if  no  BLINT  model,  next  line  is  commented  out.  If  no  yield  strength  correction 

*  factor  ‘corr’  is  omitted 

* 

blint  1  soft  1  hard  2  wsl  0.084455  wbl  0.084455  fnc=0.0  corr 
mix  3 

endep 

♦ 

tracer 
block  1 

add  0. 0.  to  0.  -10.1345  n  10 

endtracer 

* 

edit 
block  1 
noexpanded 
endblock 

endedit 

* 
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eos 

* 


*  NOTE:  EOS  properties  from  cth  mgrun  library 

* 

*  NOTE:  5083-H131  Aluminum  eos  approximated  with  6061-t6  Aluminum 

*  density  reduced  to  reflect  that  for  5083  Aluminum  from 

*  Metal’s  Handbook. 

* 

matl  mgrun  eos=6061-t6_al  r0=2.66  cs=0.534e6  s-=1.4 
gO-1.97  cv-1.07ell 

mat2mgruneos-tungsten_nir0-18.16  cs-0.403e6  s-1.237 
g0=1.67  cv=1.66el0 

endeos 

* 

*eor*  cthin 
* 

nose  shape  tests:  ogive  v=1286  target-5083  al  finite 
* 

control 
tstop  300.e-6 
rdumpf  3600 
cpshift  999. 

endcontrol 

* 

restart 

*  file=‘rsctr 
time-O.e-6 

endr 

IN 

cellthermo 

mmp 

ntbad-99999 

endc 

* 

convct 

♦ 

*  NOTE:  if  KE  Sim  convection-0,  if  MV  Sim  convection-1 

* 

convection-0 

interface-high_resolution 

endconvct 

* 

edit 

shortt 

time  0.  dtfrequency  150.e-6 
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ends 

longt 

time  0.  dtfrequency  600.e-6 
endl 
plott 

time  0.  dtfrequency  50.e-6 
endp 
histt 

time  0.  dtfrequency  0.3e-6 
htracerl 
htracer2 
htracer3 
htracer4 
htracerS 
hti:acer6 
htracer? 
htracerS 
htracer9 
htracerlO 
endh 

ende 

* 

boundary 
bhydro 
block  1 
bxbot  0 
bxtop  2 
bybot  2 
bytop  2 
endb 
endh 

endb 

* 

fracts 
pressure 
pfracl  -4.5e9 
pfrac2  -35.0e9 
pfrnix  -1.0e20 
pfvoid  -1.0e20 

endf 

* 

*eor*  pltinp 
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Intentionally  left  blank. 
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Appendix  G: 

Table  of  Simulation  Results 
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Intentionally  left  blank. 
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n  - II 

II - -  KE-NB-NC-P  Siim  vs.  5083  Aluminum  Targets  || 

65  g  Hemi-Nose  Rods  (kl-5) 

II  63  g  Ogival-Nose  Rods  (kl-4)  II 

Striking  Vebdty 
(tafs) 

Residual  Velocity 

Comments 

Striking  Velocity 
(m/s) 

Residual  Velocity 
(m/s) 

Comments 

1147 

726 

eroded,  exoestivefy 
dinned 

1600 

1375 

eroded 

1283 

930 

eroded,  exoesaivet' 
dinned 

- 

65  K  Hemi-Nose  Rods  CkMS) 

1  63  g  Ogival-Nose  Rods  (kl-5)  If 

Striking  Velocity 
(m/s) 

Residual  Velocity 
- - (m/s) 

Comments 

Striking  Velocity 
- (B^)  

Residual  Velocity 
(m/s) 

Comments  | 

1147 

927 

rigid 

1286 

1043 

rigid  1 

1201 

1041 

eroded 

1399 

1186 

rigid  1 

1283 

1021 

eroded 

1534 

1339 

rigid  1 

1600 

1439 

II -  KE-BL-NC-P  Sims  vs.  i 

083  Aluminum  Targets 

^  1 

65  e  Hemi-Nose  Rods  (id-7) 

63  g  Ogival-Nose  Rods  (id-6)  I 

Striking  Velocity 
(m/s) 

Residual  Vek>city 
- (art) - 

Comments 

Striking  Velocity 
(m/s) 

Resklual  Velocity 
(m/s) 

Comments 

1147 

828 

eroded,  excessively 
deformed 

1286 

1123 

nose  deformed  slight^ 

1201 

870 

eroded,  excessively 
deformed 

1399 

1292 

nose  deformed 

1283 

1007 

eroded,  excessive^ 
deformed 

1534 

1326 

eroded 

1600 

1420 

eroded 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  DEFENSE  TECHNICAL 
INFORMATION  CENTER 
DTICDDA 

8725  JOHN  J  KINGMAN  RD 
STE0944 

FT  BELVOIR  VA  22060-6218 

1  HQDA 

DAMOFDQ 
DENNIS  SCHMIDT 
400  ARMY  PENTAGON 
WASHINGTON  DC  20310-0460 

1  DPTY  ASSIST  SCY  FOR  R&T 
SARDTT  F  MILTON 
RM  3EA79  THE  PENTAGON 
WASHINGTON  DC  20310-0103 

1  OSD 

OUSD(A&T)/ODDDR&E(R) 

JLUPO 

THE  PENTAGON 
WASHINGTON  DC  20301-7100 

1  CECOM 

SP  &  TRRSTRL  COMMCTN  DIV 
AMSELRDSTMCM 
H  SOICHER 

FT  MONMOUTH  NJ  07703-5203 

1  PRINDPTYFORTCHNLGYHQ 
USARMYMATCOM 
AMCDCGT 
MFISETTE 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  PRIN  DPTY  FOR  ACQUSTN  HQS 
USARMYMATCOM 
AMCDCGA 
D  ADAMS 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  DPTY  CG  FOR  RDE  HQS 
USARMYMATCOM 
AMCRD 

BG  BEAUCHAMP 

5001  EISENHOWER  AVE 

ALEXANDRIA  VA  22333-0001 


1  INST  FOR  ADVNCD  TCHNLGY 

THE  UNIV  OF  TEXAS  AT  AUSTIN 
PO  BOX  202797 
AUSTIN  TX  78720-2797 

1  USAASA 

MOASAI  WPARRON 
9325  GUNSTON  RD  STE  N3 19 
FT  BELVOIR  VA  22060-5582 

1  CECOM 

PM  GPS  COLS  YOUNG 
FT  MONMOUTH  NJ  07703 

1  GPS  JOINT  PROG  OFC  DIR 
COL J CLAY 

2435  VELA  WAY  STE  1613 
LOS  ANGELES  AFB  CA  90245-5500 

1  ELECTRONIC  SYS  DIV  DIR 
CECOM  RDEC 
JNIEMELA 

FT  MONMOUTH  NJ  07703 

3  DARPA 
L  STOTTS 
JPENNELLA 
BKASPAR 
3701 N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

1  USAFSMC/CED 
DMA/JPO 
MISON 

2435  VELA  WAY  STE  1613 
LOS  ANGELES  AFB  CA 
90245-5500 

1  US  MILITARY  ACADEMY 

MATH  SCI  cm  OF  EXCELLENCE 
DEPT  OF  MATHEMATICAL  SCI 
MDN  A  MAJ  DON  ENGEN 
THAYER  HALL 
WEST  POINT  NY  10996-1786 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRLCS  ALTP 
2800  POWDER  MILL  RD 
ADELPHI MD  20783-1145 


85 


NO.  OF 

£Q£1£S  ORGANIZATION 


1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRLCSALTA 
2800  POWDER  MILL  RD 
ADELPHI MD  20783-1145 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRLQLL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 


ABERDEEN  PROVING  GROUND 

4  DIRUSARL 

AMSRLaLP(305) 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

rOPTRS  ORGANIZATION 


6  CDR  US  ARMY  ARDEC 
AMSTAARAEEWW 
E BAKER 
CCHIN 
RFONG 
J  PEARSON 
J  WALSH 
TECH  LIB 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  CDR  US  ARMY  ARDEC 

AMSTA  AR  AET  M  TECH  LIB 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  CDR  US  ARMY  ARDEC 
AMSTA  ARFS 
EANDRICOPOULOS 
PICATINNY  ARSENAL  NJ 
07806-5000 

3  CDR  US  ARMY  MERDEC 
AMSMERDSTWF 
L CRAFT 
D LOVELACE 
MSCHEXNAYDER 
REDSTONE  ARSENAL  AL 
35898-5250 


1  DIR  US  ARO  WASH 

AMXRO  W  K  A  BANNISTER 
RM8N31 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 


3  DIR  US  ARO 

J CHANDRA 
KIYER 
TECHUB 

)  PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

A 

3  DIR  US  ARMY  DARPA 

J  RICHARDSON 
TECH  INFO 
B  WILCOX 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 


4  CDR  US  ARMY  COE 
JBALSARA 
T  BLEVINS 
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